Growing Surfaces with Anomalous Diffusion - 
Results for the Fractal Kardar-Parisi-Zhang Equation 
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In this paper I study a model for a growing surface in the presence of anomalous diffusion, also 
known as the Fractal Kardar-Parisi-Zhang equation (FKPZ). This equation includes a fractional 
Laplacian that accounts for the possibility that surface transport is caused by a hopping mechanism 
of a Levy flight. It is shown that for a specific choice of parameters of the FKPZ equation, the 
equation can be solved exactly in one dimension, so that all the critical exponents, which describe 
the surface that grows under FKPZ, can be derived for that case. Afterwards, the Self-Consistent 
Expansion (SCE) is used to predict the critical exponents for the FKPZ model for any choice of 
the parameters and any spatial dimension. It is then verified that the results obtained using SCE 
recover the exact result in one dimension. At the end a simple picture for the behavior of the Fractal 
KPZ equation is suggested, and the upper critical dimension of this model is discussed. 



The Kardar-Parisi-Zhang (KPZ) equation |l| for sur- 
face growth under ballistic deposition was introduced as 
an extension of the Edwards- Wilkinson theory Q- The 
interest in the KPZ equation exceeds far beyond the in- 
terest in evolving surfaces because of the following rea- 
sons: (a) The KPZ system is known to be equivalent to a 
number of very different physical systems. Examples are 
the directed polymer in a random medium, Schrodinger 
equation (in imaginary time) for a particle in the pres- 
ence of a potential that is random in space and time and 
the important Burgers equation from hydrodynamics • 
(b) The second reason, that is more important, to my 
mind, is that it serves as a relatively simple prototype of 
non-linear stochastic field equations that are so common 
in condensed matter physics. 

The equation for the height of the surface at the point 
r and time t, h (r, t) , is given by 



^ = vV 2 h + ±(Vhf+ V (r,t), 



(1) 



where automatically the constant deposition rate is re- 
moved and r\ (r, t) is a noise term such that 



(V(r,t))=0 

(r, (r, t) rj (f, t)> = 2D 6 (f - f*)5(t- t') 



(2) 



As can be seen in eq. Q the basic relaxation mecha- 
nism in the KPZ equation is a Laplacian term that results 
from nearest neighbor hopping in the growing surface. 
During the last years, there has been a g rowing interest 
in other relaxational mechanisms |22 | -|25 | . namely subd- 
iffusive diffusion that seems to appear in the context of 
charge transport in amorphous semiconductors |26l |27| , 
NMR diffusometry in disordered materials |28l. and the 
dynamics of a bead in polymer networks 29] . The the- 
oretical effort to account for such phenomena led to the 



formulation of the celebrated fractional Fokker-Planck 
equation (FFPE) |22|. This equation includes a frac- 
tional (Riemann-Liouville) operator instead of the stan- 
dard derivative of the Fokker-Planck equation. Recently, 
Mann and Woyczynski |j| have suggested that in order 
to account for experimental data, namely experiments in 
which impurities were present on the growing surface Q , 
a modification of the KPZ equation has to be considered. 
They used the observation that the presence of an impu- 
rity can act as a strong trap for an adatom migrating at 
room temperature, to conjecture that this process corre- 
sponds to Levy flights between trap sites. This conjec- 
ture then served as a justification for the introduction of a 
fractional Laplacian into the continuum equation of the 
growing surface as another relaxation mechanism. Ac- 
tually, the fractional Laplacian dominates the standard 
Laplacian in the KPZ equation in the scaling regime (i.e. 
in the large scale limit), so that the standard Laplacian 
can be ignored from the beginning. To summarize this 
exposition and to be more specific, the equation they 
eventually suggest to describe the growing surface in the 
presence of self-similar hopping surface diffusion is the 
fractal KPZ (FKPZ) equation given by 



f = -(A^) + ^(VM 2 

W 2 



V (r, t) , 



(3) 



where A Q = — (—A) ' (or in Fourier space (A Q T) (q) = 
, 2 W 2 



(wry 



\q\ J T (q)) is the fractional Laplacian, and in our 

context it is more convenient to choose a = 2 — p with 
< p < 2 (where the special case p = corresponds 
to the standard KPZ equation). In addition, ij (r,t) is a 
noise characterized by 



(T) (?,t)ri(?,t)) = 2D \r-f'\ 2a ' d 5{t-t l ) 



(4) 
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Actually Mann and Woyczynski Q discussed the specific 
case of white noise that corresponds to a = in the 
last equation, but since the more general case does not 
require special efforts we discuss the FKPZ problem with 
spatially correlated noise. Furthermore, an exact solution 



2 



is possible only for a special case with correlated noise, 
so that the more general discussion is also interesting on 
that basis. 

The proposed FKPZ equation generalizes the FFPE 
equation mentioned above in that it is a field equation, 
rather than an equation for a single degree of freedom. 
Therefore, it is understood that such a generalization is 
essential in order to account for the dynamics of a whole 
medium experiencing anomalous diffusion, and not just 
an artificial problem. Obviously, the technical mathe- 
matical difficulties to be overcome in this nonlinear case 
are formidable in comparison to those for the linear frac- 
tal kinetic equations. 

However, in their paper Mann and Woyczynski Q were 
not able to predict the critical exponents that describe 
the surface that grows under FKPZ. But before I make 
any new statements about this model let me briefly sum- 
marize the various quantities of interest. 

A very important quantity of interest is the roughness 
exponent a that characterizes the surface in steady state. 
The roughness exponent is usually defined using W (L, t) 
the roughness of the surface (that is defined as the RMS 
of the height function h (r , t) , in a system of linear size 
L). Then, in terms of W (L, t), a is given by 

W (L, t) oc L a . (5) 

Another important quantity of interest is the growth 
exponent (3 that describes the short time behavior of the 
roughness W (L, t) (with flat initial conditions) 

W(L,t) cx t 13 . (6) 

Finally, I introduce the dynamic exponent z that de- 
scribes the typical relaxation time scale of the system 
(i.e. the dependence of the equilibration time on the size 
of the system) 

t x oc L z . (7) 

It is well known [2| that these three exponents are not 
independent and that under very general considerations 
one should expect the following scaling relation 

z = a/(3 (8) 

(this relation is a direct consequence of the Family- Vicsek 
scaling relation [f|. 

In addition, for the KPZ equation there is another scal- 
ing relation, that comes from a symmetry of the equation 
under infinitesimal tilting of the surface (this symmetry is 
just the famous Galilean invariance of the Burgers equa- 
tion) 

a + z = 2. (9) 

It can easily be checked that this symmetry holds in the 
case of FKPZ as well, because the fractal Burgers equa- 
tion is evidently invariant under Galilean transformation 
(see ref. 0, eq. (7.8)). Therefore, the last scaling rela- 
tion is also relevant in the this discussion. Hence both 



scaling relations, reduce the number of unknown expo- 
nents to just one (out of the three we started with). In 
some cases an extra scaling relation is possible (for ex- 
ample, in the case of the KPZ equation with long-range 
noise - see refs. till 113); so that the exponents can be 
obtained exactly by power-counting. Employing the ter- 
minology of the Dynamic Renormalization Group (DRG) 
approach this can be explained by saying that certain 
terms in the dynamic action do not renormalize, and so 
an extra scaling relation arises (see |3q). This kind of 
solutions are naturally available also in the present prob- 
lem. Whenever the exponents are obtained due to such 
an extra condition it will be specifically pointed out. 

In this paper I show that for a specific choice of 
parameters of the FKPZ equation (namely p = 2a), 
the equation can be solved exactly in one dimension, 
so that all the critical exponents can be derived eas- 
ily for that case. Afterwards, in order to give a more 
complete picture (i.e. for any dimension d, and any 
spatial correlation index a) I apply a method devel- 
oped by Schwartz and Edwards @-@ (also known as 
the Self-Consistent-Expansion (SCE) approach). This 
method has been previously applied successfully to the 
KPZ equation. The method gained much credit by be- 
ing able to give a sensible prediction for the KPZ criti- 
cal exponents in the strong coupling phase, where many 
Renormalization-Group (RG) approaches failed, as well 
as Dynamic Renormalization Group (DRG) 0,13 (actu- 
ally, it can be shown that the strong coupling regime is 
inaccessible by DRG even when it is used to all orders 
H3> El)- ^ is then verified that the results obtained 
using SCE recover the exact result in one dimension. 

As mentioned above, for the specific case p = 2a in 
one dimension an exact solution can be found for the 
FKPZ problem using the Fokker-Planck equation associ- 
ated with its Langevin form (i.e. eq. ©)• This particu- 
lar choice of parameters corresponds to a situation where 
the fractional exponent p equals the exponent 2a that 
describes the decay of spatial correlations in the noise. 
Since kind this exact solutions is familiar in the KPZ 
community I will simply state the final results given by 
a = 1/2 and z = 3/2 (using the scaling relation ©). 
Notice that these critical exponents extend the classical 
one dimensional KPZ exponents for non zero p's and cr's, 
as the classical KPZ case corresponds to p = 2a = 0. 

Now, the Self-Consistent Expansion (SCE) is applied 
in order to learn about the behavior of this system in 
more general contexts (other dimensions and cases where 
p ^ 2a). SCE's starting point is the Fokker-Planck form 
of FKPZ, from which it constructs a self-consistent ex- 
pansion of the distribution of the field concerned. 

The expansion is formulated in terms of <p q and uj q , 
where <f> q is the two-point function in momentum space, 
defined by 4> q = (h q h- g ) s , (the subscript S denotes 
steady state averaging), and uo q is the characteristic fre- 
quency associated with h q . It is expected that for small 
enough q, <p q and uj q are power laws in q, 



3 



Aq 



and 



Bq z 



(10) 



where z is just the dynamic exponent, and the exponent 
r is related to the roughness exponent a by 



r-d 



(ii) 



The main idea is to write the Fokker-Planck equation 
dP/dt = OP in the form dP/dt = [O + O x + 2 ] P, 
where Oq is to be considered zero order in some 
parameterA, Oi is first order and Oi is second order. The 
evolution operator Oq is chosen to have a simple form 



Oo 



E 



()lln 



D. 



d 

1 dh- a 



qhq 



where D q /uj q 



Note that at present 4> q and ui q are not known. I obtain 
next an equation for the two-point function. The expan- 
sion has the form (j) q — 4> q + c q {{4> P } , because the 
lowest order in the expansion already yields the unknown 
(f) q . In the same way an expansion for oj q is also obtained 
in the form uj q = uj q + d q {{4> P } , {w p }). Now, the two- 
point function and the characteristic frequency are thus 
determined by the two coupled equations 



First, there are two kinds of weak-coupling solutions 

- both with a dynamic exponent z = 2 — p (they are 
called weak-coupling because they are exactly the solu- 
tions obtained in the case of the Fractal Edwards- Wilkin- 
son equation, see |3j). Now, when the spatial correlations 
of the noise are relevant (i.e., when a > 0) I obtain the 
solution r = 2 — p + 2(7- provided that d > 2 — 3p + 2a. 
But if the spatial correlations of the noise are not relevant 
(i.e., when a < 0) I obtain the simpler solution T = 2 — p 

- provided that d > 2 — 3p. 

The second type of solutions is strong coupling so- 
lutions that obey the well known scaling relation z = 
d+ 2~ T obtained from eq. I|14|) (this scaling relation is 
just the above-mentioned scaling relation a + z = 2 
that is naturally obeyed by our analysis). The first 
strong coupling solution is determined by the combina- 
tion of the scaling relation and the transcendental equa- 
tion F (r, z) — 0, where F is given by 



({^ p } , K|) = o 



and 



d, H}) = 0(12) 



These equations can be solved exactly in the asymp- 
totic limit to yield the required scaling exponents gov- 
erning the steady state behavior and the time evolution. 
Working to second order in the expansion, one gets the 
two coupled integral equations 



m 2 \ ^ ly - L q£m lvJ - ImqH'mY' i 
9 q U)q + UJ£ + UJ m f-* 1 UJ„ + UJ£ + Ul, 



M q e m Mi mq (j) m (j) q 



M q l m M m £ q (j)£(f) q 

LO„+U>e+ 0J m 



0. 



(13) 



and 



2E M < 



qtm 



Ml mq 4> rn + M m £ q (j)e 
UJ e +UJ m + UJ q 



0, (14) 



where D q — D q~ 2cr , v q — vq 2 ~ p and M q £ m = 

^= (i- rfij S q j +m . In addition, in deriving eq. (|14f> I 

have used the Herring consistency equation 01 • I n f ac t 
Herring's definition of ui q is one of many possibilities, 
each leading to a different consistency equation. But it 
can be shown, as previously done in pj, that this does 
not affect the exponents (universality). 

A detailed solution of equations l|13fl and (|14fl in the 
limit of small q's (i.e. large scales) in the line of refs. Q- 
@ yields a rich family of solutions that I shall describe 
immediately. 



F(T,z) 



, t ■ (e - t) 

d d t , v ' x 

t z + \e - tj +1 

(e-t)r T + e- (e-t) |e-i]~ 

d d t _\tfiz3_ t -T \ g _ tT r ,(i 5 ) 
t* + \e-i\ +1 1 1 V ^ 

and e is a unit vector in an arbitrary direction. 
This solution is valid as long as the solutions of the 
last equations satisfy the following condition T > 
max{ d+4 3 + 4g ,^,d + 2 / o}. 

It turns out that for d = 1 the equation F (T, z (T)) = 
is exactly solvable, and yields T = 2 and z = | (it can 
be checked immediately by direct substitution). In this 
case the validity condition reads 2 > max { 5 3 CT , 1 + 2p} 
or cquivalcntly < p < 1/2 and a < 1/4. By using eq. 
(flip I translate the results into a = \ and z = | that 
are precisely the exact result presented above. 

It should also be mentioned that for d > 2 such an 
exact solution in closed form cannot be found, and one 
has to solve numerically the equation F(T, z (T)) = 0. 
For convenience I denote the numerical value of this so- 
lution by Tq (d). For example, in two dimensions I obtain 
r (2) = 2.59. 

The second strong coupling solution is obtained by 
power-counting, and it is relevant when a > (i.e. 
when the spatial correlations of the noise are relevant), 
d < 2 + 2cr and F (T, z (T)) < (the last condition turns 
out to be equivalent to the condition r > Tq (d) because 
from this value of T and on the function F (T, z (T)) is 
negative). This power counting solution can be writ- 
ten in closed form and is given by z = d+4 3 ~ 2cr and 

r = d+4 3 + 4 ' T = z + 2cr. 

The third strong coupling solution (that is in some 
sense the only "genuine" FKPZ solution, in the sense 
that it is the only solution that is dramatically influenced 
by the fractional Laplacian, and at the same time is not 
a solution of the Fractal Edwards- Wilkinson equation) is 
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also a obtained by power counting. More specifically, it 
is determined by the combination of the scaling relation 
z = d+ 2~ T an d the extra relation d+4—2Y—z = 2—p—T. 
This solution can be written in closed form as Y = d + 2p 
and z = 2 — p. It turns out that this phase is relevant 
when d > max {2 — 3p, 2 — 3p + 2a}. In addition, it is 
needed that F (r, z (Y)) is positive. Therefore, this solu- 
tion is possible only when Y < Tq (d) . 

The following table summarizes all the possible phases 
found in this paper 



z 


r 


Validity 


2-p 


2 - p + 2a 


a > and d > 2 - 3p + 2a 


2-p 


2-p 


a < and d > 2 - 3p 


2-p 


d + 2p 


< p < 1 and 


d+4-2tr 
3 


d+A+Acr 
3 


d > max {2 - 3p, 2 - 3p + 2a} 
a > 0, d < 2 + 2a and 


d+4-r (d) 

2 


r (d) 


d + 4 + 4cr > 3r (d) 

To (d) > max { d+4 + 4CT , d + 4 , d + 2p} 



TABLE I: A complete description of all the possible phases 
of the FKPZ problem, for any value of d, p and a. The first 
two columns give the scaling exponents z and F for a partic- 
ular phase, and the third column states each phase's validity 
condition. Note that To (d) is the numerical solution of the 
transcendental equation F (T, z) = with the scaling relation 
z = d+4 F T - if such a solution exists. 

In order to gain more insight on this system, it might 
be interesting to specialize to two extreme cases: namely 
p = and a ^ vs. p ^ and er = 0. The first case 
(namely p = and a ^ 0) corresponds to the local KPZ 
problem with long-range noise. This problem has been 
studied in the past using various methods - for example: 
DRG HmUl, Mode-Coupling and SCE 0. All 
methods agree on the basic picture that for a big enough 
noise exponent (er) one obtains a power-counting strong- 
coupling solution, given by z = d + i ^ 2cr (z [ s the dynamic 
exponent). The controversy between the different meth- 
ods is over the values of the scaling exponents for smaller 
values of er, and on the critical value errj that separates be- 
tween the two phases. Not surprisingly, the results given 
here agree with the previous SCE result presented in ref. 
□ 

The second case (namely p ^ and a = 0) corresponds 
to the Fractal KPZ problem with white noise that is the 
original equation suggested by Mann and Woyczynski Q . 
Specializing the general picture of phases presented in 
Table I above to the case a = yields new phase di- 
agrams, namely a separate phase diagram for every di- 
mension d. For example, in two dimensions there are 
three possible phases: I. The standard KPZ phase given 
by T (2) = 2.59 that is valid for < p < 0.295. II. 
The weak-coupling phase given by Y = 2 — p that is valid 
for any p > 0. III. The "third" strong coupling solution 
given by Y = d + 2p that is possible for < p < 1. In 
the first phase the dynamic exponent is z — 1.705, while 
the other two phases share the same dynamic exponent 



of z = 2 — p. The possible phases in two dimensions are 
presented in Fig. 1. 
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FIG. 1: The values of the scaling exponents (a) z and (b) V 
as a function of the fractional parameter p, for uncorrelated 
noise (a = 0) in two-dimensions. 

The described two-dimensional (and white-noise) sce- 
nario is quite typical and appears in other dimensions 
as well - namely there are usually three possible phases 
that are possible for different values of a with possible 
phase transitions between them (as in the usual KPZ sce- 
nario - the phase transition is controlled by the strength 
of the dimensionless coupling constant). More precisely, 
this picture extends up to the upper critical dimension of 
the original KPZ problem {d^ c pz ) where the first phase 
(phase I in Fig. 1) disappears (see the discussion below). 

At this point, the full picture of possible solutions 
might seem too complex, so I want to suggest the fol- 
lowing simple interpretation for the behavior of the Frac- 
tal KPZ equation. If you remember, the starting point 
of this discussion was the introduction of the fractional 
Laplacian into the KPZ equation. This immediately im- 
plies that faster relaxations are now possible (faster in 
the sense that a smaller dynamical exponent is expected) 
when compared with the Edwards- Wilkinson equation. 
However, it is well known that already the KPZ nonlin- 
earity introduces faster relaxations (at least for dimen- 
sions lower than the upper critical dimension). There- 
fore, in the FKPZ system, the dynamics is controlled 
by the fastest " component" : if the dynamical exponent 
of the classical KPZ system is smaller than 2 — p then it 
dominates, otherwise the new fractional dynamical expo- 
nent z = 2 — p controls the dynamics. However, at this 
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point the picture gets a little bit more complicated (just 
like in the classical KPZ case), namely there are several 
possible phases with this new fractional dynamical expo- 
nent z = 2 — p, and the transition between these phases 
is controlled by the strength of dimensionless coupling 
constant. 

These conclusions have an important implication for 
the upper critical dimension of the FKPZ model (i.e. 
the dimension above which the dynamical exponent is 
the same as that of the linear theory). Namely, it is 
turns out that the FKPZ equation always has an up- 
per critical dimension that is determined by the relation 
r o {dic FZ ) = d lc PZ + 2 P ( or alternatively, using the 
roughness exponent a - the upper critical dimension is 
the dimension where qkpz (4 Z ) = P)- Notice, that 
this result does not dependent on the ongoing debate 
over the existence of the upper critical dimension for the 



KPZ system (see ^2" EH)- Actually, it merely requires 
that the roughness exponent of the classical KPZ system 
becomes arbitrarily small in higher dimensions - an as- 
sumption that is generally accepted (only Tu [2l| had a 
different conjecture). 

As one can see, the results obtained using the Self- 
Consistent Expansion (SCE) are quite general, and cover 
all possible values of the relevant parameters (p and a) 
as well as dimensions. It is also easily verified that these 
results recover the exact result obtained at the beginning 
of this paper for the case p = 2a. This situation, suggest 
that the SCE method is generally appropriate when 
dealing with non-linear continuum equations. 
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